BNL-69354 
FERMIL AB-Conf-02/227-T 
IFT-P.070/02 

1 



Possible Origin of RHIC R ou t/R sld HBT Results 

(N ; 

§ ■ Sandra S. Padula^ 
(N 



Oh- 

OS 



On 
O 



o 

=5 



X 



a Instituto de Ffsica Teorica, Universidade Estadual Paulista, Sao Paulo - SP, Brazil. 

The effects of opacity of the nuclei together with a blackbody type of emission along the 
system history are considered as a means to explain the ratio R ou t/R S id observed by STAR 
and PHENIX collaborations at RHIC. Within our model, no flow is required to explain 
the data trend of this ratio for large surface emissivities. 
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The unexpected results presented by STAR[[IJ in the previous Quark Matter meeting on 
7r ± 7r ± HBT, later confirmed by PHENIX[ Q], regarding the decrease of the ratio R ou tl ' Rsid 



C\| ■ for increasing Kt, has been challenging explanations since that time. Hydrodynamic 
models and microscopic-based simulations usually predicted the opposite behavior for 
increasing pair momentum. 

Motivated by this challenge we proposed a simple model (see [ [| for details) to try 
and understand the unforseen decrease of the outwards radius (along the direction of the 
average transverse momentum of the pair of pions, Kt) relative to the sidewards one (i.e., 
orthogonal to Kt)- Two were the main ingredients of this attempt. The first was to 
consider the particle emission of a blackbody type, radiating from the external surface of 



the system during its entire evolution. The second was to treat this system as an opaque 
source. The other basic ingredients were not unusual. The system produced in a heavy ion 
collision is supposed to be formed above the critical temperature, T c , at the time r , in a 
Quark-Gluon Plasma (QGP) phase. Its temperature gradually decreases while it expands 
and for the sake of simplicity, the expansion is considered to be only along the longitudinal 
direction. After this initial stage, lasting about (r c — r ), where r c correspond to the on-set 
of the phase-transition at the temperature T c , the mixed phase begins, during which the 
temperature remains constant with time. The mixed phase continues for a longer period, 
ending after an elapsed interval (r^ — r c ). Then, the system converted into a gas of pions 
(no resonances are included) expands further, until the decoupling temperature, Tf, is 
reached. At this point, the system is quite dilute, since most of the particles have already 
been evaporated from its surface. Thus, at the time 77, the system is supposed to decouple 
in an instantaneous volumetric emission. No complex mechanism for the hadronization of 
quarks and gluons is considered in detail at this point, although hadronization must take 
place. In other word, in first approximation, the evaporation of "gluons" and "quarks" 
(as hadronized pions) from the external surface of the system is considered in the same 
way as emission of pions, except for the number of degrees of freedom. 

♦Theoretical Physics, FERMIL AB, P.O. Box 500-MS106, Batavia, IL 60510, USA. 
^Nuclear Theory Group, Brookhaven National Laboratory, Upton, NY 11973, USA. 
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The Bjorken hydrodynamical model[ || is assumed to describe the system during its 
entire evolution, i.e., since it is formed and until it breaks up. This is supplemented with 
a blackbody type of radiation from the surface of the matter from its formation. The 
emitting source is supposed to be opaque, in a generalized version of opacity proposed by 
Heiselberg and Vischerf || at CERN energies, later followed by Heinz and Tomasik [ [|. 

To compute the emitted spectrum and the two particle distribution function, the Co- 
variant Current Ensemble formalism, [|7j, is adopted. In this formalism, the two particle 
correlation function can be written as C(k u k 2 ) = P ^)p k Jk 2 ) = 1 + G(^M)G{k 2 MV where 
Piiki) and P 2 (ki, k 2 ) are, respectively, the single particle distribution and the probability 
for simultaneous observation of two particles with momenta k\ and k 2 - The average and 
the relative momentum of the pair are defined as K = (k x + k 2 )/2 and q = kx — k 2 . 

The emitted energy as well as the total entropy associated to each stage of the system 
evolution can be estimated. Just for a brief illustration, I write down the emitted energy 
as a function of time in the initial stage by considering the emission by an expanding 
cylinder of transverse radius Rt and length h, in the time interval r and r + dr, as 
dE in = — koT^Iii R^hdr — ^aT^nR^dh, where the first term comes from the blackbody 
type of energy radiated from the external surface of the cylinder, and the second term 
results from the mechanical work due to its expansion. The k factor was introduced to 
take into account that the system has some opacity to surface emission. The constant a is 
the Stefan-Boltzmann constant and is proportional to the number of degrees of freedom 
in the system. By integrating this expression the energy density can be obtained as 
6i„ = eo(^-) 3 e r t^ t . From this expression it can be seen that the multiplicative factor, 

e R r y , appears in addition to that coming from the Bjorken picture. The variation 
of the temperature in the initial stage, i.e., prior to the beginning of the phase transition, 
follows immediately as T(r) = T (^)^e 2R r^ T T °\ The instant corresponding to the 
beginning of the mixed phase, r c , when T c = 175 MeV is reached, that one corresponding 
to its end, at Th, as well as and the decoupling time, 77, at T c = 150 MeV, are given by 

Finally, the initial values of the temperature, T and the formation time, To, are related 
to the initial entropy, So and to the input number of particles, J\f (chosen to match the 
average experimental pion multiplicity per unit of rapidity at RHIC, M ~ 1000), by 
S = TAf = [(g 9 + g q ) x (|)gT 3 ] rcR 2 T r ; r ~ ^(GeV- 1 ) ~ °-g(fm), where T = 3.6, 
as estimated by the entropy per particle (S w /N n ) of a pion gas at freeze-out. Then, 
T ~ 411 MeV and r ~ 0.160 fm. The degeneracy factors, g g + g q , account for the gluon 
and quark (antiquark) degrees of freedom (g g + g q = 37 for two quark flavors). In the case 
of pions, the degeneracy factor is g n = 3. I illustrate in the table below the time variables 
for two different assumptions on the emissivity, k. I also write the estimated fraction of 
the particles emitted from the surface during the period r < r < Tf, S/J\f to t, relative to 
the total number of produced particles, Ntot, as well as the remnant portion at freeze-out, 
V/Aftot, then emitted from the entire volume. 
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In order to check how the spectra estimated within our model behave as compared to 
data (PHENIX minimum bias [ ||), I plot its predictions on the single-inclusive distri- 
bution in Fig. 1(a). The estimates and discussions presented here are restricted to the 
central rapidity region, i.e., yi = (which implies that ki L = 0, and, consequently, Kl = 
and qi = 0). In Fig. 1(b) I show the results for the ratio R out /R S id vs. Kt, together with 
the preliminary STAR[ [IJ (filled triangles) and PHENIX[ Q (filled circles) data for both 
7r + 7r + and 7r~7r~ interferometry. 
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Figure 1. (a): The prediction based on our model for the transverse momentum distribu- 
tion of emitted pions is shown. The points are from the minimum-bias data from PHENIX 
Collaboration. The curves correspond to emissivity k = 0.5 and to k = 1, without the 
inclusion of any transverse flow. Both cases describe data on spectrum well in the low pion 
momentum region, up to about k ir « 1 GeV/c. (b): The results for the ratio R ou t/R S id of 
the outwards by the sidewards radii are shown within our model. The ratio corresponding 
to full emissivity (k = 1) agrees very well with data within the experimental error bars 
(shown in the plot), whereas the 50% emissivity case is completely excluded by data. The 
parameters used are given in the text and in the table, corresponding also to T = 411 
MeV, T c = 175 MeV, Tf = 150 MeV, and the transverse radius, ~ 7 fm/c. 

As seen from the plots in Figure 1, our results were extremely successful in describing 
both sets of data for k = 1, but the curve corresponding to k — 0.5 is away above the data 
limits, suggesting that there should be high emissivity along the system history in order 
to explain the data trend. The model also describes the typical source radii reasonably 
well, but not the Kt dependence of these radii (see [ |J for details). This suggests that 
the time variation of the emitting radius and the introduction of transverse flow may play 
a significant role, [ [|-[ |i~H| . If there is a time variation of the various radii, this will be 
correlated with the typical momentum scale of emitted particles, since the earlier is the 
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time, the hotter are the particles. Our model is sensitive to such variation since emission 
from the hot surface is allowed at early times. Also, a proper treatment of the decoupling 
is not included in our computations, which would affect the results, although it might 



also suggest modification in the treatment of decoupling, [ pl]|,[ |12|. 

The principal reason why such a small ratio of R ou t/ Rside vs. is obtained within 
our model is probably due to a combination of two effects. The first is that the surface is 
opaque, and whatever is emitted from the surface will have a small value of this radius. 
The second effect is that black body radiation by partons is allowed when the surface is 
very hot. This allows a much larger contribution from surface emission than is typical of 
what happens in hydrodynamical simulations, where particles are emitted by Cooper-Frye 
decoupling from a surface at very low temperature. In fact about 80% of the emission 
comes from the surface in our model. The fact that so many particle are emitted from the 
surface at early times also means that the longitudinal decoupling time in this computation 
is significantly shorter than would be the case for hydrodynamic simulations. I should 
add that many of the features of the model proposed and discussed here are embodied in 
the hydrodynamic computations of Heinz and Kolb[ The main difference lies in the 
treatment of an essential ingredient, the emissivity of the surface. 
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